Momentum dependence of the symmetry potential and nuclear reactions induced by 

neutron-rich nuclei at RIA 
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Effects of the momentum-dependence of the symmetry potential in nuclear reactions induced by 
neutron-rich nuclei at RIA energies are studied using an isospin- and momentum-dependent trans- 
port model. It is found that symmetry potentials with and without the momentum-dependence but 
corresponding to the same density-dependent symmetry energy Esym{p) lead to significantly differ- 
ent predictions on several i7sHm (p)-sensitive experimental observables. The momentum-dependence 
of the symmetry potential is thus critically important for investigating accurately the equation of 
state (EOS) and novel properties of dense neutron-rich matter at RIA. 
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The rapid advance in technologies to accelerate ra- 
dioactive beams has opened up several new frontiers 
in nuclear sciences 0, Particularly, the high en- 

ergy radioactive beams to be available at the planned 
Rare Isotope Accelerator (RIA) and the new accelera- 
tor facility at GSI provide a great opportunity to ex- 
plore the EOS and novel properties of dense neutron- 
rich matter[l|l|10|. The energy per particle E{p, S) 
in asymmetric nuclear matter of density p and relative 
neutron excess S = (p„ — Pp)/{pp + Pn) is usually ex- 
pressed as E{p,5) = E{p,S = 0) -f E^yn,{pW + 0{5^), 



cause only very recently was the momentum-dependence 
of the symmetry potential given in a form practically 
usable in transport model calculations [i^. In this work, 
we study for the very first time effects of the momentum- 
dependence of the symmetry potential within an isospin- 
and momentum-dependent transport model for nuclear 
reactions induced by neutron-rich nuclei at RIA ener- 
gies. It is found that symmetry potentials with and 
without the momentum-dependence but corresponding 
to the same symmetry energy Esym{p) lead to signifi- 
cantly different predictions on several sensitive probes of 



where Esym^p) is the density-dependent nuclear sym- the Esym{p). Moreover, these observables are more sen 



^sym 

metry energy. The latter is among the most impor- 
tant and yet very poorly known properties of dense 
neutron-rich matter For instance, it is important 
for Type II supernova explosions, for neutron-star merg- 
ers, and for the stability of neutron stars. It also de- 
termines the proton fraction and electron chemical po- 
tential in neutron stars at (3 equilibrium. These quan- 
tities consequently determine the cooling rate and neu- 
trino emission flux of protoncutron stars and the possi- 
bility of kaon condensation in dense stellar matter 0, ■ 
In nuclear reactions induced by neutron-rich nuclei, the 
Egymip) reveals itself through dynamical effects of the 
corresponding symmetry potentials acting differently on 
neutrons and protons. Based on isospin-dependent trans- 
port model calculations, several experimental observables 
have been identified as promising probes of the Esym if) , 
such as, the neutron/proton ratio the neutron- 

proton differential flow 0, 0| , the neutron-proton cor- 
relation function Ti*! and the isobaric yield ratios of light 
clusters [l5j. However, in all existing transport models 
the momentum-dependence of the symmetry (isovector) 
potential has never been taken into account. Effects of 
the momentum-dependence of the symmetry potential in 
nuclear reactions are thus completely unknown, although 
effects of the momentum-dependence of the isoscalar po- 
tential are well-known ^Id. .IT^ .IQj . This is mainly be- 



sitive to the variation of Esym if) in calculations with the 
momentum-dependent symmetry potentials. 

That the momentum-dependence of the single particle 
potential is different for neutrons and protons in asym- 
metric nuclear matter is well-known and has been a sub- 
ject of intensive research based on various many- b ody 
theories using non-local interactions, see e.g., ref. |2l| 
for a review. Guided by a Hartree-Fock calculation us- 
ing the Gogny effective interaction, the single nucleon 
potential was recently parameterized as |20| 



U{p,6,p,T) = A^t^ + AA 
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In the above t — 1/2 (—1/2) for neutrons (protons) and 
r 7^ r'; cr = 4/3; frir^p) is the phase space distribu- 
tion function at coordinate r and momentum p. The 
parameters AutAi,B,Ct,ttCt,t' and A listed in ref.jioll 
were obtained by fitting the momentum-dependence of 
the U{p, S,p,t) predicted by the Gogny Hartree-Fock 
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and/or the Brueckner-Hartree-Fock calculations, satura- 
tion properties of symmetric nuclear matter and the sym- 
metry energy of 30 MeV at normal nuclear matter den- 
sity po — 0.16/ f TV?. The compressibility of symmetric 
nuclear matter Kq is set to be 211 MeV. The momentum- 
dependence of the symmetry potential steams from the 
different interaction strength parameters Cr,T' and Cr^r 
for a nucleon of isospin r interacting, respectively, with 
unlike and like nucleons in the background fields. More 
specifically, Cuniike = -103.4 MeV while Cuke = -11.7 
MeV. 

The parameter x was introduced to reflect the largely 
uncertain behavior of Esy^ip), specifically its poten- 
tial part E^y^{p). The latter corresponding to a; = 1 
(denoted by MDI(l) which gives the same Egy^ip) as 
the default Gogny interaction) and x — Q (MDI(O)) 
together with the kinetic contribution E^y^{p) — 
?i^/67Ti„(37r^/9/2)^/'^ are shown in Fig. 1. The variation 





60 








40 






> 


20 































E 
S 

in 


-20 




^«^""' - 












-40 








-60 







P'Po 

FIG. 1: The density dependence of the potential and kinetic 
parts of nuclear symmetry energy. 

of Esym{p) resulted by changing from the MDI(l) to 
the MDI(O) parameter set is well within the uncertain 
ran ge p redicted by various many-body theories, see e.g., 
refs.|22L E3 . li^ . The potential contribution to the sym- 
metry energy can be parameterized as 

E^ymip) = 3.08 + 39.6u- 29.2^2 -k5.68u^-0.523m^ (2) 

for the MDI(l) (Gogny) and 

E^ymip) = 1-27 + 25.4u - 9.31u2 + 2.171*^ - 0.21^" (3) 

for the MDI(O), where u = p/po is the reduced nucleon 
density. 

We implemented the single particle potential in eq. 
in an isospin-dependent transport modelQ. Since the 
C-terms in the single particle potential depend insepa- 
rably on the density, momentum and isospin, to inves- 
tigate effects of the momentum-dependence of the sym- 
metry potential we shall compare results obtained using 



eq. n] with those obtained using a single nucleon po- 
tential U„omsip,S,p,T) = Uo{p,p) + Usym{p,5,T) that 

has almost the same momentum-dependent isoscalar part 
U[){p,'p) as that embedded in eq. and a momentum- 
independent symmetry potential Ugym {p, 5, t) that gives 
the same Esym{p) as eq. ^ The Usym{p,S,T) can be ob- 
tained from Usym{p,5,T) = dWsym/dpr, where Wgym is 
the isospin-dependent part of the potential energy den- 
sity Wsym = EPy^{p) ■ p-S^. Using the parameterizations 
for Efy^{p) in eqs. El and 13 we obtain 

Usynr'Hp, (5, t) = 4r(5(3.08 -I- 39. 6u - 29. 2^^ + 5.68^^ 
-0.52u4) - ^2(3. 08 + 29.2m2 - 11.4^3 + 1.57^4) (4) 

for the MDI(l) parameter set and 

U^^^^'°ip, 6, t) = 4r(5(1.27 + 25.4u - 9.31^^ + 2.17^^ 
-0.21u4) - (5^(1.27 + 9.31^2 - 4.33^^ + 0.63u'') (5) 

for the MDI(O) parameter set. To identify reliably effects 
of the momentum-dependence of the symmetry poten- 
tial without much interference from density effects the 
Unomsip,S,P,T) should lead to almost the same reaction 
dynamics and evolution of density profiles as the single 
particle potential in eq. ^ Both of them are mainly de- 
termined by the isoscalar potential for which we select 
the original MDYI interaction^^! 

Uo{p,p} = -110.44m + 140.9wi-2'* 

fd^W IMl (6) 

PO J ^ l + (p-p')V(l-58p°^)^'^ ^ 

where is the Fermi momentum. The compressibility 
Kq for this interaction is 215 MeV. We compared numer- 
ically Uo{p,p) in eq. Elwith {Un{p,S,p) + Up{p,S,p))/2 
obtained from eq. 1. They are indeed very close and both 
agree with the nucleon optical potential data. We also 
verified numerically that the potential UnomsiPTS,P,T) 
constructed this way leads to about the same reaction 
dynamics and the evolution of density profiles as the po- 
tential in eq. ^ 

The strengths of the symmetry potentials with and 
without the momentum-dependence but corresponding 
to the same Egymip) are compared in Fig. 2 by exam- 
ining the difference between neutron and proton poten- 
tials ([/„ — Up)/2. The symmetry potential without the 
momentum-dependence is higher in magnitude and has 
generally steeper slopes than the momentum-dependent 
one for p/po < 2.3 with the MDI(l) parameter set and 
at all densities for the MDI(O) parameter set. Moreover, 
the strength of momentum-dependent symmetry poten- 
tials decreases with the increasing momentum. Thus 
the difference between the symmetry potentials with and 
without the momentum-dependence is larger for nucle- 
ons with higher momenta. Several experimental observ- 
ables are known to be sensitive only to the symmetry 
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FIG. 2: .Strengths of the symmetry potentials with and with- 
out (sohd) the momentum-dependence as a function of density 
as measured by the difference between neutron and proton 
potentials. 



potential but not to the isoscalar potential. These ob- 
servables are mainly neutron-proton differential or rel- 
ative quantities^, where effects of the isoscalar po- 
tential with or without the momentum-dependence are 
largely canceled out. These observables are also in- 
sensitive to the in-medium nucleon-nucleon (NN) cross 
sections 0, 0, Q- We use here the isospin-dependent 
experimental NN cross sections. Assuming the symme- 
try potential is momentum-independent, these observ- 
ables were previously proposed as promising probes of the 
density-dependence of the symmetry energy. We com- 
pare in the following several such observables calculated 
with the U{p,5,p,t) in eq. ^ and the Unomsip,S,P,T) 
corresponding to the MDI(l) parameter set. Conclusions 
based on calculations using the MDI(O) parameter set are 
the samejl^. 

Shown in the upper window of Fig. 3 is the average 
isospin asymmetry Sfreeiu) of free nucleons as a function 
of rapidity y from 1400 events of ^^'^Sn -|-i24 reactions 
at a beam energy of 400 MeV/nucleon and an impact 
parameter of 5 fm. We define here free nucleons as those 
with local baryon densities less than po/8. The value of 
Sfreeiy) reflects mainly the degree of isospin fractiona- 
tion between the free nucleons and the bound ones at 
freeze-out. It is also influenced slightly by the produc- 
tion of more tt~ than tt'^ mesons in the reaction. The 
momentum-independent symmetry potential leads to sig- 
nificantly higher value of Sfreeiy) than the momentum- 
dependent one. Moreover, the difference tends to increase 



with rapidity. At midrapidity the predicted Sfreeiv) val- 
ues are close to the value expected when a complete 
isospin equilibrium is established among all target and 
projectile nucleons. These features are what one expects 
from the strength of the symmetry potentials as shown 
in the upper window of Fig. 2. The generally repulsive 
(attractive) symmetry potential for neutrons (protons) 
around po causes more neutrons (protons) to be free 
(bound) . The momentum-independent symmetry poten- 
tial is higher and steeper than the momentum-dependent 
one and the difference between them increases with mo- 
mentum, it thus leads to the higher Sfreeiy) values espe- 
cially at higher rapidities. 
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FIG. 3: Isospin-asymmetry (upper window) and neutron- 
proton differential flow (lower window) of free nucleons as 
a function of rapidity. The solid (dashed) lines are calculated 
with the momentum-independent (-dependent) symmetry po- 
tential. 



Shown in the lower window of Fig. 3 is the neutron- 
proton differential flow F^_p(y) = Y^f^.f^ip^Wi) / Niy), 
where Wi = 1(— 1) for neutrons (protons) and Niy) is 
the total number of free nucleons at rapidity y. The 
differential flow combines constructively effects of the 
symmetry potential on the isospin fractionation and the 
collective flow. It has the advantage of maximizing ef- 
fects of the symmetry potential while minimizing ef- 
fects of the isoscalar potential|0. Compared to the 
momentum-dependent symmetry potential embedded in 
eq. ^ the momentum-independent symmetry potential 
U^ym^^{p,S,T) makes not only more neutrons free but 
also gives them higher transverse momenta in the reac- 
tion plane because of its higher magnitude and steeper 
density slopes. As a result, the differential flow with 
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the momentum-independent symmetry potential is sig- 
nificantly higher. 

Complementary information about how the symmetry 
potential depends on the momentum can be obtained by 
studying the ratio of neutrons to protons as a function 
of their transverse momentum pt- Shown in Fig. 4 is this 
ratio around the midrapidity within \ycms/ybeam\ < 0.3. 
The overall rise of the ratio at low pt is due to the 
Coulomb force which shifts protons from low to higher 
energies. It is seen that the difference between the pre- 
dicted ratios increases with pt, reflecting the increas- 
ingly larger difference between the symmetry potentials 
with and without the momentum-dependence for nu- 
cleons with higher momenta as shown in Fig. 2. The 
high Pt nucleons are thus more useful for studying the 
momentum-dependence of the symmetry potential. 
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FIG. 4: Tfie ratio of free neutron to proton multiplicity as a 
function of transverse momentum at midrapidity. The solid 
(dashed) line is calculated with the momentum-independent 
(-dependent) symmetry potential. 



Are the observables studied above still sensitive to 
the variation of Egym (p) when the momentum-dependent 
symmetry potential is used? To answer this question 
we have compared calculations using the MDI(l) and 
MDI(O) parameter sets with the momentum-dependent 
and -independent symmetry potentials. Using the po- 
tential in eq. these observables are actually more sen- 
sitive to the variation of Esym{p) than the calculations 
with the momentum- independent symmetry potentials. 
For instance, the slope of the differential flow F^_p{y) 
at y=0 changes from -1.9 to 21.4 MeV/c using eq. 
by changing from the MDI(l) to the MDI(O) parame- 
ter set, while it changes from 26 to 39 MeV/c with the 
momentum-independent symmetry potentials. The net 
change due to the variation of the Egym ifi) is thus larger 
with the momentum-dependent symmetry potentials. 

In conclusion, the momentum-dependence of the sym- 



metry potential is found to play an important role in 
heavy-ion collisions induced by neutron-rich nuclei at 
RIA energies. Symmetry potentials with and without the 
momentum-dependence but corresponding to the same 
symmetry energy Egymip) lead to significantly differ- 
ent predictions on several experimental observables that 
were previously identified as promising probes of the 
Egymip)- With the momentum-dependent symmetry po- 
tential these observables are more sensitive to the change 
of Esym{p)- The momentum-dependence of the symme- 
try potential is thus critical for investigating accurately 
the EOS of dense neutron-rich matter at RIA energies. 
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